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Introduction and Summary 

There is a large numbers of solvable models in statistical mechanics where the bulk 
free energy and possibly other quantities such as order parameters and the correlation 
length can be calculated exactly. Mostly these models are two-dimensional, and only a 
few three-dimensional examples are known. The first such model was the two-dimensional 
Ising model, solved by Onsager in 1944 [1]. Then, it took some two decades to realize 
that two important ingredients of the Ising model — the star-triangle relation (STR) and 
the resulting commutativity of the transfer matrices — can be used to solve the other 
two-dimensional models [2,3]. 

There is now quite rich theory of solvable two-dimensional lattice models. There are 
several different (but related) methods, mostly based on the commutativity of the transfer 
matrices, and thereby on the Yang-Baxter relation which have been developed (for a review 
see [4,5]). 

Can we find three-dimensional models with commuting transfer matrices? It is known 
that the tetrahedron relation [6] replaces the Yang-Baxter relation as a commutativity 
condition [7,8] for the three-dimensional cubic lattice. This relation contains thousands 
of distinct algebraic equations, and, obviously, it is very difficult to solve them. The only 
non-trivial solution obtained by direct analysis of these equations known so far is that of 
Zamolodchikov [6,9]. In fact, even in the simplest two state spin case one has to solve 2^^ 
equations, instead of 2^ in two dimensions. The symmetry properties could slightly reduce 
this number, but the jump in complexity is still enormous. Therefore one would like to 
find an alternative approach to the commutativity which would be based on more simple 
algebraic relations. 

In the present paper we discuss one such scheme. We consider the solvable interaction- 
round-a-cube model [10] on the cubic lattice with A^-valued spins (A^ > 2) at each site. To 
within a minor modification of the boundary condition it is equivalent to the s/(n)-chiral 
Potts model [11] (this model has an interesting history which can be traced in [12-17]). 
On the other hand it can be regarded as a multistate generalization of the Zamolodchikov 
model [6], reducing to it when N = 2. The Boltzmann weight function of eight corner 
spins around a cube have a very special form, such that by introducing an auxiliary center 
spin for each cube the model can be viewed as an Ising-type one on the body centered 
cubic lattice with two- and three-spin interactions only. The corresponding Boltzmann 
weights are of course not arbitrary, and obey a number of simple relations. The most 
complicated among them is a "restricted star-triangle relation" (RSTR) which relates 
two- and three-spin weights. Remarkably enough this relation is a simplified version of 
the two-dimensional star-triangle relation (for the two-state spin case it can be obtained 
simply by specialization of rapidities in the STR of the 2d Ising model). 

Starting from these relations we derive symmetry properties of a Boltzmann weight 
function of the elementary cube of the lattice under transformations from the symmetry 
group of the cube. As in the Zamolodchikov model this Boltzmann weight function depends 
on three parameters which can conveniently be chosen to be the dihedral angles 61,62, ^3 
between three "rapidity planes" passing through the cube. The symmetry properties are 
entirely consistent with such geometric interpretation of these angles. 
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The implication of the symmetry properties is twofold. First, they contain the three- 
dimensional star-star relation [10,18]^ which results in the commutativity of the transfer 
matrices. The other symmetry relations imply certain symmetry properties of the partition 
function per site for the lattice infinite in the all three directions. These two features — 
commutativity and symmetry — allow us to find the partition function per site extending 
the N = 2 Zamolodchikov case calculations of [18] to the general N case. The result is 
quite surprising: 

log«Ar(^l,^2,^3) = ^^^^^log«2(^l,^2,^3) 

where kat is the partition function per site for the A^-state model. In addition we show 
that there exists a three-dimensional free boson (Gaussian) model whose partition function 
per site is given by 

log«B(^l,^2,^3) = -21og«2(^i,^2,^3) 

This means that kn can be expressed as a (rational) power of some free boson determinant. 
It would be very interesting to understand the reason for this fact. 

The sequence of our working is illustrated on Fig.l. We hope this will assist the reader. 

1. Formulation of the model 

1.1 The interaction-round-a-cube model 

Consider a simple cubic lattice £ of M sites with the periodical boundary conditions 
in each direction. At each site of C place a spin variable s taking N > 2 distinct values 
s = 0, . . . , — 1, and allow all possible interactions of the spins within each elementary 
cube. The partition function reads 

^=X1 n ^^i^o\si,S2,S3\si,S2,Ss\so) (1.1) 
spins cubes 

where sq, ■ ■ ■ ,so are the eight spins of the cube arranged as in Fig. 2, and V{sq S2, 
^sl-Si, S2, S3 1 So) is the Boltzmann weight of the spin configuration sq, . . . , sq. Note that 
the spins sq, si, S2, S3, si, S2, S3, sq correspond respectively to a, 6, c, d, e, /, g, h in 
Fig. 2 of reference [10]. The product is over all elementary cubes in C 

Taking the lattice to have m horizontal layers and letting 0j denote all spins in layer 
i, one can rewrite (1.1) as 

^ ^■■■^T^i'l>2T<j>2H---T<i>rn<l>i=T^T'^ , (1-2) 

where T is a layer-to-layer transfer matrix whose elements are the products of all the V 
functions of cubes between two adjacent layers. Clearly, T depends on the Boltzmann 
weight function so we can write it as T{y). 

^ Interestingly, this relation corresponds to the "central inversion" of the cubic lattice. 
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Obviously, we cannot solve the model (1.1) for an arbitrary Boltzmann weight function 
V. What we apparently can do is to solve a restricted class of models whose transfer 
matrices form commuting families, i.e., such that any two transfer matrices T{V) and 
T{y') belonging to the same family commute, 

[T{V),T{V')] = . (1.3) 

In all known cases this commutativity condition can always be reformulated as a 
relation for the local Boltzmann weight functions V and V . In two dimensions the corre- 
sponding relation is well known: it is the Yang-Baxter equation for the local Boltzmann 
weights. This equation involves three different Boltzmann weight functions. A straight- 
forward generalization [7,8] of this construction for a 3-dimensional lattice leads to the 
tetrahedron equation [6] , involving four different Boltzmann weight functions. Even in the 
simplest two- valued spin case this equation involves thousands of distinct relations and, 
obviously, is very diflficult to solve. The only known non-trivial solution obtained from the 
direct analysis of these relations is the Zamolodchikov's one [6,9]. 

In this paper we present a somewhat different approach to the commutativity in three 
dimensions, which could, however, be applied to the Zamolodchikov model as well. We con- 
sider the interaction-round-a-cube model of ref [10] where the Boltzmann weight function 
of eight corner spins around a cube, V, has a very special form, such that by introducing 
an auxiliary center spin for each cube the model can be viewed as an Ising-type one on the 
body centered cubic lattice with two- and three-spin interactions only. The correspond- 
ing 2-spin and 3-spin Boltzmann weights obey a number of simple relations including "a 
restricted star-triangle relation" (RSTR) which allow us to prove the commutativity of 
the transfer matrices and calculate the partition function per site by use of the symmetry 
properties of the model. 

Before ending this subsection note that there are simple transformations of V that do 
not change (1.1). For instance, if we multiply V{so,...,s) by F{so, S2, si, s^)/ 
F{'^i, ss, So, S2) then each horizontal face of C acquires an F-factor from the cube below 
it, and a canceling (1/F) factor from the cube above. The effect on the transfer matrix 
T is to apply a diagonal similarity transformation. Provided F is the same for T{V) and 
T(V') the commutation relation is unaffected. This is an example of a "face-factor" trans- 
formation. Similarly, one could apply "edge-factor" and "site-factor" transformations that 
leave (1.1) and (1.3) unchanged. 

1.2 The Boltzmann weights 

Let a, b, c be integers and v be a complex variable. Introduce the following notation 
uj = exp{27ri/N) , a;^/^ = exp(m/7V) , (1.4) 

7(a, 6) = o;'^'' , $(a) = (cc;^/2)a(iv+a) ^ 

A{v) = {l-v'^y/'^ . (1.6) 
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Note that 



7(a, b) = -f{b, a) = 7(a + N, b) , 

$(a + iV) = $(a) . (1.7) 
Define a function w{v, a) sucli tliat 

^ = (Awrn(i-.^-.)- . (1.8) 

Obviously, the LHS of this equation is a multivalued function of v, but it is a singlevalued 
function of a point {v, A(t>)) on the algebraic curve 

^;^ + A^ = l . (1.9) 

Below we shall interpret w in this latter way, suppressing, however, an explicit dependence 
on the phase of A{v) in the arguments of w. 

Now, fix four complex parameters p, p' , q, q' , and define 

vi = q'/{ujp') , V2 = q'/p , v^=p/q , VA=p'/q . (1.10) 

Obviously 

LVV1V4 = V2 Vs . (1-11) 

The Boltzmann weight function of the model is given by the equations (2.9), (2.10) 
of ref.[10]. Taking into account equation (1.18) (given below) and omitting equivalence 
transformation factors (which do not change the partition function) one can write it in the 
form 

V{so\si, S2, S3\si, S2, S3\so\\vi, V2, V3, V4) = 

^ ^ y^^ w{v2, si-S2 + a)w {vs, S3 - sp - a) ^{sp, (7)7(50, a) ,^ 
^ w(vi, S3-so + a)w {v4, S1-S2- a) 7(52, (^his2, o") 

where p is a normalization factor depending on f 1, . . . , 1)4. Note that the Boltzmann weight 
function (1.12) describes a very special type of interaction of eight spins around the cube. A 
typical cube with its centre spin, cr, is shown in Figure 2. There are three- spin interactions 
on the shaded triangles, described by w{v, a) or by l/w{v, a) in (1.12). There are also 
two-spin interactions such as 7(so, cr), or 1/7(^2, cr) associated with the edges linking cr 
to So, So, S2, S2 (these edges are denoted by heavy lines in Figure 2.). In addition, the 
three-spin and two-spin interaction weights w and s obey a number of important relations 
which we now discuss. These relations, rather than the explicit form of w and 7 given by 
(1.5) and (1.8), enable us to execute the program stated in the introduction and calculate 
the partition function per site in the thermodynamic limit. 
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1.3 Restricted star-triangle relation. 

Consider the following automorphism of the curve (1.9) 

{v, A{v)) ^ (v, Aiv)) , (1.13) 

— , A(v) = —. (1.14) 

UJV V 

One can easily show that 7(0, b), $(a) and w{v, a) satisfy the following properties 

7(a, 6 + c) = 7(a, b) 7(0, c) (1.15) 

N-l 

J2 l{a,b)-f{-b,c) = N5ac (1.16) 
6=0 

$(a + 6) = $(a) $(6) 7(a, 6) (1.17) 

w{v, a) tv{v, —a) 



w{v, 0) !(;({;, 0) 



= $-^(a) (1.18) 



where f is related to v by equation (1.14). Moreover these functions satisfy two less trivial 
relations. For the moment let vi, V2, V3, denote four arbitrary complex variables, then 

I _ w{v2,a-e) w{v'2,-b) w{v2/iujvi), a) 
Fi{vi, V2\a, b)= y — = ifiivi, V2) ^, , 1.19 



F2 ^3, ^4 a, h)=Y, ^ = ^3, ^4 , , , / r , 1.20 

^ 'u;(v4, a-£) w\y'^, -b)w{v^lV2„ a) 

where and <^2 are scalar functions and 

^ UVxA{v2) A{V2) 

, A(-^i) , V2A{y2l{^vx)) , . 

"^=A(^ ' ^^"^^ = A(.,) ' 

, t'4A(t)3) , A('i;4/'i;3) 



= 



V3A(^;4) A(^;4) 



^(4)= , (l-21c) 
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Each of the relations (1.19), (1.20) is a corollary of another and the properties (1.15) - 

(1.18) . The phases of A(w4/z;3) and /S.{v2/ iojvi)) can be chosen arbitrarily since they 
cancel out the RHSs of (1.19) and (1.20). 

The reader may have noticed that relation (1.19) (or (1.20)) is a particular case of 
the usual star-triangle relation. Indeed, the LHS of (1.19) is the sum of the product of the 
three functions (each depending on two spins) while the RHS is the product of three such 
function. Unlike the usual star-triangle relation there is some asymmetry in the LHS of 

(1.19) : the function 7 does not depend on any continuous parameters. It is quite possible 
that (1.19) is a particular case of a more general relation and 7 is just a limiting value 
of a more complex function. In fact, this is exactly so ior N = 2 when (1.19) and (1.20) 
can be obtained be a specialization of rapidities in the the star-triangle relation of the 
two-dimensional Ising model [1]. This is the reason why we call the (1.19), (1.20) as the 
"restricted star-triangle relations" . 



2. The symmetry relations 



2.1 The cube symmetry 



Consider the cube C shown in Fig. 2. The eight spins sq, . . . , -^o at the corners of C 
can be grouped into four ordered pairs 

dj = {sj,Sj) , j = 0, ...,3 (2.1) 

corresponding to four diagonals do, . . . ,dj, which we assume to be oriented and directed 
from Sj to Sj, J = 0, . . . , 3. The spatial symmetry group of the cube Q{C), (consisting of 
all possible reflection and rotations which map the cube to itself) has a structure of the 
direct product 

g{C) = 02x84. (2.2) 

of the cyclic group of order 2, C2, generated by the central inversion, P (which reverses 
the directions of all four diagonals do, ... , ds) and the symmetric group, S4, of order 24, 
consisting of the transformations which permute the diagonals preserving their directions.^ 
Let 

S0123 = (sq S2, S3I Si, S2, S3I So) (2.3) 

denote the sequence of the spins sq, . . . ,so, corresponding to their basic arrangement at 
the corners of C as shown in Fig. 2. 

Obviously, the transformations from ^(C) induce some permutations of the spins 
So, . . . , So- According to the above discussion of the structure (2.2) of g{C) these transfor- 
mations map the sequence (2.3) either to Sijki or to s^h, 

Sijke = {si \sj, Sk, Si] Sj, Sfc, se\ s) (2.4) 
Sijke = {si \sj, Sfc, s^l Sj, Sfc, se\ s) (2.5) 



^ For the even permutation the corresponding transformation is a pure rotation of the 
cube, while for the odd permutation it is a rotation followed by the central inversion, P. 
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where {i, j, k, £) is some permutation of (0, 1, 2, 3). 

Consider two elements of G{C), specifying their action on the spins. Let {i, j, k, £) be 
any permutation of (0, 1, 2, 3), and let R e G{C), 

Rsijke = sekij 1 Rsijke = sekij , R^ = \ (2.6) 

denote the 90° rotation around the axes passing through the centers of the top and bottom 
faces of C, while T e g{C), 

Tsijk£ = Sikje , Tsijk£ = Sikj£ , = 1 , (2.7) 

denotes the reflection with respect to the plane passing through the corners of C occupied 
by the spins sq, S3, Sq, S3 in Fig. 2. 

One could easily check that these two elements R and T generate the whole group 
G{C). In particular, the central inversion, P, 

Psiju = Sijkl 5 Psijki = Sijki , = 1 , (2.8) 
can be expressed as 

P = {RTf . (2.9) 

2.2 The angle parameterization 

Let ^1, 92-I Os denote angles of a spherical triangle and ai, 02, 03 denote three sides of 
this triangle opposite to the angles 61, $2, ^3. Define the related variables 



CKo = (6*1 + 6*2 + 6*3 - 7r)/2 , cij = 6'i-Q;o , 

Pq = (27r - ai - 02 - 03)72 , Pi = TT - - ai , 



(2.10) 



for i = 1, 2, 3. Choose 9i, 62, Os, so that 6*1, 6*2, 6*3, ao, . . . ,a2, Po, . . . , P3 are all real, 
between and tt. Further, define (taking real positive values of roots) 



(2.11) 



= [sin(^,/2)]'/^ , a = [cos(^,/2)]'/^ , 

T, = [tan(^i/2)]'/'^ , Zi = eMm/N) , 
for z = 1, 2, 3 and 

Ui = eMiPi/N) , Q = [cos(V2)]'/2 , z = 0,...,3 . (2.12) 
Now parameterize p, p', q, q' in (1.10) as follows^ 

p=z^^T^ , = a;-i/%-iTr^ , q = T^^ , q' = u-^'^T2, 



^ This parameterization has been obtained by a formal generalization of eqs.(4.19),(4.24) 
of ref.[10] for arbitrary values of N. 



then^ 

V3 = z~^TiT2 , V4 = u~^/'^z^^T2/Ti . 
Now choose the phases of A(f 1), . . . , A{v4) such that 

Aiv^)=Ss/{CiC2Us) , A{v2) = CsVi/{SiC2) 

AK) = S3Us/{C^C2) , A(^4) = Cs/{S^C2Ui) 



(2.13) 



(2.14) 



Aiv^/vs) = A{v2/{u;ui)) = S^^ (2.15) 

With these definitions the Boltzmann weight function (1.12) can be regarded as a 
function of the three independent variables 6*1, 6*2, 6*3 or, equivalently, of the four dependent 
variables cio, . . . , 03 constrained by the relation 

CKo + CKi + Q!2 + Q!3 = TT (2.16) 

so we can write the LHS of (21) as V{so, . . . , so||q;o, • • • > c^s)- 

Note, that the angles ^1,^2,6*3 can be viewed as the dihedral angles between three 
"rapidity planes" passing through the cube exactly as it is in the Zamolodchikov model. 



2.3 The normalization of the weights 

Now we like to fix the normalization factors w{v, 0) in (1.8) and p in (1.12). 
First let us choose w{v, 0) such that 

N-l 

Yl w{v, a) = 1 . (2.17) 

a=0 

With this normalization define 

D+{v) = {detN\\w{v,a-b)\\y/'' , 
D.{v)^{detN\\l/w{v,a-b)\\y/'' . 

Also, set 



(2.18) 



5+ = (detAr||7(a, 6)11)1/^ , 
<S_ = (detAr||l/7(a, 6)il)'/^ 



(2.19) 



^ For N = 2 the variables vi, . . . , here differ from those given by (7.19) of [18] merely 
by negating V3 and V4. 
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Regarding a, b in (1.19), (1-20) as matrix indices running the values 0, . . . , A^— 1, and taking 
the determinants of the both sides of these equations one gets by using (2.17)-(2.19) 

(2.20) 

¥^2(^3, ^4) = D_(v4)S+/D+(v'^) , 
where v'l, v'^ are given by (1.21). Explicit calculations with the equations (1.8) give 

D±iv) = c± {v/A{v))^^-^^/'^ 
S+S_ = N 



(2.21) 



where c± are (inessential) constants. Note, in particular, that when fi , . . . , f 4 are param- 
eterized by (2. 13)- (2. 15) we have 

■ a \ {N-i)/N 
sm 6'2 ^ 



MV1,V2)MV3,V4) = N (2.22) 

Vsm Os J 

Further, set the normalization factor in (1.12) as 

p = (202(7V-i)/iV/^ ^ (2.23) 
1 \ 1/2 



2^=\^-sm0sj /(C0C1C2C3) (2.24) 

where cq, ■ ■ ■ ,cs are given by (2.12). 

We want to calculate the free energy, or equivalently the partition function per site 

K = Z^l^ . (2.25) 

Note that when ckq = ci;2 = we have from (2.10) - (2.13) 

= = ^^3 = ^^4 = , 2^ = 1 , (2.26) 

while from (1.8) 

t(;(0, a) = t(;(0, 0) , Va . (2.27) 
Substituting (2.26), (2.27) into (1.12) and using (1.15), (1.16) one gets 

V = bs,--s,^s.--s, . (2.28) 

Ignoring irrelevant boundary contributions, it follows that for M large 

K = 1 , when ckq = 0:2 = . (2.29) 
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2.4 The symmetry properties of the Boltzmann weights 



The definition (1.12) of the Boltzmann weight function V is obviously rather asym- 
metric with respect to the orientation of the elementary cube (see Fig. 3). We shall see, 
however, that despite this visual asymmetry, the Boltzmann weight function (1.12) has a 
remarkable hidden symmetry. Below we shall show that up to the equivalence transfor- 
mation factors, which do not affect the partition function, the Boltzmann weight function 
remains unchanged upon the permutations the corner spins sq, . . . ,so induced by sym- 
metry transformations of the cube complemented by corresponding transformation of the 
variables vi,V2,vz, V4 (or, conveniently, the variables ckq, . . . , 0:3). Obviously, it is enough 
to prove this just for two transformations (2.6) and (2.7) since they generate the whole 
symmetry group of the cube. 

From (1.3), (1.18), (2.17) it follows that 

w{v, 0)w{v, 0) = e^'^^^'-i)/^ (2.30) 

is a constant independent of v. Let us of interchange ai with 0:2 leaving ao,as intact. 
From (2.13) this is equivalent to the replacement of vi , V2,V3, by vi , V4, V3,V2 respectively. 
Using (2.30) and noting that (2.23) remains unchanged one can easily trace the effect of 
this interchange on the weight function (1.12) 

l^(soi23||ao,ai,a2,Q;3) = ^7^^ — =4v^(so2i3| l^o, "2, "i, "3) (2-31) 

^{Sl - S2) 

where we have used the short notations (2.4) for the spin configuration of the corner spins. 
This gives the transformation law of V under the reflection T, (2.7). 
Further, using (1.15), (1-16) rewrite (1-12) in the form 



s{so- S3,So- S2)p ( S{ss + Ss, H) 



X 



s{so-S3,so-S2)N ^ is(so + so,/i) (2 32) 

Fi{vi,V2\so + S1-S2- S3,S2 - So + At)-F2 (t'3, ^'4|so + Si - S2 - S3, S2 - Sq - At)| 

Applying now (1.19), (1.20) and taking into account (2.14), (2.15), (2.20)-(2.24) one obtains 
V{soi23\\ao, ai, a2, as) = 

S(-S() - .So - :S2) H'(('4/('3, .S() + .Si - /?2 - ^3) . , , x (2.33) 

-7= = T— ^ -, _ , _ T 1/ (S32OI "3, "2, "0, ai) 

s(So - S3, So - S2) W{V4/V3, So + Si - S2 - S3) 

which gives the required transformation of V under the rotation, R, (2.6). Combining now 
eqs.(2.32) and (2.33) and using (1.18), (2.30) for the functions w in (2.33) one gets the 
relation corresponding to the i?T-transformation 

l^(soi23||Q!o,ai,Q:2,Q;3) = 

^{so)^{ss) s{s2,So- Ss)w{v2V3,Si + Ss-So-S2)^,,_ . (2-34) 

^7=-T^7=-T-7 I r V^(s3io2 q;3, ai, "o, "2) 

$(So)$(S3) S(S2, So - S3) W{V2V3, Si S3 - So - S2) 
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where the argument of tu's can be written as 

V2V3 = q'/q = u;-'^^Ti. (2.35) 
Remembering eq.(2.9) and iterating (2.34) three times one obtains 

^(■Soi23||ao,ai,a2,Q;3) *^(so) s{s2,so - S3) s(so,S3 - S2) s(s3,S2 - So) 



(2.36) 



V{soi23\\ao,ai,a2,a3) $^(so) s(s2,So - S3) s(so, S3 - S2) s(s3,S2 - Sq) 
W{V2V3, Si + S3 - So - S2) W{V5, Sj + ^2 - Sp - S3) w{v4/v3, Sp + Si - S2 - S3) 
W{V2V3,SI + S3 - So - S2) W{V5, Si + S2 - Sq - S3) w{v4/v3,So + Si - S2 - S3) 

where 

^4A('^i)A('t;3) ,-1/2^2 A/ X ^2 .^o7N 
V3A(t>2)A(^;4) 

One can check that (2.36) is exactly the three-dimensional star-star relation conjectured 
previously (eq.(6.1) of ref.[10]). 

As we remarked before the angles 9 1,62, O3 can be interpreted as the dihedral angles 
between the three rapidity planes rigidly connected with the cube. Then, from geometric 
considerations these angles should be very simply transformed by the cube symmetry group 
Q{C). Namely, the related variables ctQ, cti, a2, 03 given by (2.10) should just permute for 
any transformation from G{C). This is entirely consistent with (2.31), (2.33). 

3. Partition function 

3.2 Factorization and commutativity. 

Consider two successive layers of C with / spins (p on the lower layer, cp' on the upper. 
In the center of each intervening cube we have the central spin a as in Fig. 3. Let (p" denote 
the set of all these a-spins between (p and (p" . Then because the top spins sq, si, S2, S3 in 
Fig. 3 interact only with one another and with a, and similarity for thee bottom spins, we 
can write the elements of the transfer matrix T as 



where 

TT w(w3,S3 - So - cT)7(so,a) 

t(;(t'4,Si - S2 - cr)7(S2,cr) 
cubes ("i 0\ 

y-f w{v2,Si-S2 + cr)'y{sQ,a) 

The products are over the all the I cubes between the two layers; for each cube so, si, S2, S3, 
si,S2,S3,so are the eight spins shown in Fig 3. 
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Obviously we can regard X^^^n as the element of a matrix X. From (3.2) this matrix 
depends on and V4, so it can be written as X{v^, v^). With similar conventions for Y , 
(3.1) implies 

T ^ p'X{v^,va)Y{vi,V2). (3.3) 

Thus p~^T factors into a product of two matrices, one dependent on ^1, ^2, ^3 only via 
and f4, the other via vi and V2- 

In [10] it was shown that transfer matrices T form two-parameter commuting families. 
More precisely, when parameterizing PiP'-iQiq' through the angles (2.13), the eqs. (1.5), 
(2.11), (2.17), (2.22) of ref. [10] imply that two transfer matrices (with different values of 
1,02,03 ) commute provided they have the same value of Oi 

[T{Oi, O2, Os), T{Oi, 0'^, 0'^)] = 0, y02, Os, 0'^, 0', (3.4) 

The proof was based on the Yang-Baxter equation for the ,s/(?i)-chiral Potts model [11, 
17,19]. Alternatively we can establish the commutativity property (3.4) directly from the 
3-dimensional star-star relation (2.36) not referring to those results. This will be done 
below in this section. 

If Oi is known then V2 can be determined from vi and from vs by corollaries of 
(2.14). 

/ OA 

V2 = o;^/^ I cot Vi, 

^ ^ (3.5) 



1/2/ (^iV^"" 

V4 = UJ I 1 cot y ) 



The text step is to show that the factorization property (3.3) remains true when the 
matrices T, X, Y are appropriately diagonalized. To ensure this, it is necessary to slightly 
modify the model. 

There are two sorts of vertical faces in C: those whose perpendiculars run in front-to-back 
direction (such as S0S2S3S1 and S3S1S0S2 in Fig. 3), and those whose perpendiculars run 
right-to-left. Call the former type "FB", the latter "RL". At the center of each FB face 
place a spin fx, with values 0, . . . , — 1. Let the spins on the front and back faces in Fig. 3 
be // and respectively. Choose them so that 

a = pi- n' (mod N). (3.6) 

Do this for all cubes in £ . If a' is the spin behind cr, and a" is the spin behind that, etc.; 
then on using the cyclic boundary conditions we observe that 

(j + (j' + a" + --- = {pi-pi') + {pi' - p") + {p" - p'") + • • • = 0. (mod N) (3.7) 

(Each |U-spin occurs twice with the opposite signs. If C has n layers perpendicular to the 
front-to-back direction then there are n u-spins on the LHS of (3.6)). Thus we can use 
(3.6) only if the sum of each horizontal front-to-back line of u-spins is constrained to be 
zero. This is merely a change of boundary conditions, and in the limit of n large it should 
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have no effect on the partition function per site k. We shall refer to the model subject to 
these constraints as the "modified model" . 

For the modified model (3.1) and (3.2) formally remains the same, but (p" is now the 
set of all /x-spins and a in (3.2) is now given by (3.6). Note also that in this case the 
matrices X, Y, T are unchanged under overall shifts of all //-spins, or all s-spins, on any 
front-to-back line. Therefore X, Y and T have non-zero entries only in the diagonal block 
with respect to the subspace invariant under all such shifts. 

For the moment, let us ignore the constraints (3.5), and work with the original model. 
The elements of the transfer matrix XY in (3.1)-(3.3) are obtained by taking the product 
over all cubes in a layer of the function F(soi23 || ao, ■ ■ ■ ■, "s) given by (1.12), (2.3). Let us 
replace the function V by V = V(soi23 || «o, • • • , cts) with S0123 given by (2.5) and denote 
the corresponding transfer matrix as T. Then with the same arguments which led to (3.3) 
one can show that 

T = p'X{vi,V2)Y{vs,V4), (3.8) 

where X, Y are similar but not identical to X, Y in (3.3). From (2.36) the functions V 
and V differ by the equivalence transformation factors only. If we take the product of 
these factors over all cubes in a layer then the result can be put in the form L(0)/L(0'), 
where (f) is the set of spins on the lower layer, 0' on the upper. 
We have therefore shown that 

[Xivi,V2)Y{v3,V4)U,^> = L{cl))[X{v3,V4)Y{vi,V2)U,4,'/L{cl)') (3.9) 

or in matrix notation 

X{VUV2)Y{V3,V4) = LX{V3,V4)Y{VUV2)L-^ (3.10) 

where L is the diagonal matrix with elements L{(f))S{(f),(f)'). It depends on di, but not on 
any other parameters. 

Now introduce the constraints (3.6). This is equivalent to introducing a factor 

n{^E-''"""'""""-4 (3-11) 

I fc=0 J 

into summand in (3.1), where a, a' ,a" , . . . are the n center spins an a horizontal front- 
to-back line, and the outer product is over l/n such lines in the horizontal layer of spins 

r- 

Let Vfe be a function that differs from V only in that an extra single spin factor cj''"' 
is put into the summand in (1.12). Inserting (3.11) into (3.1) we obtain 

L fc=o ) 

where the inner product is over n cubes in a horizontal front-to-back line and again the 
outer product is over all l/n such lines in a layer. From (1.2) 

Vkiso\si,S2,S3\si,S2,S3\so) = V{so + k\si, 82,83 + k\8i, 82, 83\so) (3.13) 
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Shifting So and S3 in (2.36) we obtain 



Vfc _ 7(fc, So - S3) W{V5, Si + S2 - Sq - S3 - k) W{V5, Si + S2- Sq- S3) V 

Vk 'y{k,so-ss)w{v5,si + S2-so-ss-k)w{v5,si + S2-so-S3)V 

where Vk = Vfc(soi23)- If we take the product of each side of (3.14) over the n cubes in 
a horizontal front-to-back line, the w and s factors cancel. It follows that we still have 
relations (3.10) in the modified model^. 

Now consider the case when V2 — = 1. From (3.5) it follows then 

vi = ^ = a;-i/2(tan ^^/N (3.15) 
UIV4 2 

while (1.12) gives 

-\ Ai5(s2 - Si,S3 - So) , . 

F(so,...,so) = ^ — = — — = ^ (3.16) 

W{V1, S3 - So - Si -I- S2)W{V4, Si - S2 - S3 - So) 

where Ai is a scalar factor. (From (2.17) Ai is infinite, but this is just a feature of the 
normalization (2.17), and can readily be removed, leaving the following arguments intact). 
The function Vk contains extra uj'^^^'^~^'-^^ factor but this cancel out of the product in (3.12), 
so in this case 

T^,^'= n V{so,...,so) (3.17) 

cubes 

where V is given by (3.16) and the product is over all cubes in a layer. Substituting (3.16) 
into (3.17) and taking into account (1.18), (2.30) one obtains 

T^,,^' ^ X2 Yl ^{s2- si,S3-so) (3.18) 

cubes 

where A2 is a constant. 

The product of Kronecker delta-functions is zero unless = 0', or if 0,0' differ only 
by overall shifts of spins on some horizontal front-to-back lines. Restricting attention to 
the subspace invariant under these shifts we then obtain from (3.3), (3.18) 

X{l,V4)Y{vi,l) ^ XI (3.19) 

where A is another (nonzero) constant, / is the unit matrix and vi and V4 are related by 
(3.15). For the modified model the matrices in (3.18) are square and hence non-singular 
and invertible^. 



^ Applying similar arguments to the product of the weight function V along one front- 
to-back line of the cube, rather that to the whole horizontal layer, one obtains precisely 
the two dimensional star-star relation of the s/(n)-chiral Potts model (eq. (3.19) of [11]). 

^ Curiously enough, it is this apparently innocuous statement that fails for the original 
model and is the reason for introducing the constraints (3.6). 
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Regard di as fixed, vi and vs as independent variables, and V2 and V4 as given by 
(3.5). Then we can suppress 9i, ^2, V4 dependence and write (3.10) as 

^(^i)^^(^3) = LX{v3)Y{vi)L-^ (3.20) 

Then (3.19) gives 

X{l)Y{x) ^ X{x)Y{l) = XI (3.21) 

where x = a;~^/^(tan^i/2)^/^ is the value of vi when V2 = 1. 

Using (3.20), (3.21) we can now perform the following transformations 



X{v3)Y{x)X{l)Yiv^) = XX{vs)Y{v^) = XL-^X{vi)Yivs)L 

L-^X{vi)Y{l)X{x)LY{vs)L = X{l)Y{vi)X{vs)Y{x) 



(3.22) 



Hence the matrix products X {l)Y{vi) and X{v-i)Y{x) commute for all values of Vi and 
v^. Assuming that these commuting matrix products can be simultaneously diagonalized 
it follows that there must exists diagonal matrices A{v3) and B{vi) and a non-singular 
matrix P (independent of the variables vi and V3) such that 

X{v3)Y{x) = PA{v3)p-\ X-^X{l)Y{vi) = PB{vi)p-^ (3.23) 

Setting Q = Y{x)P and remembering that all matrices depend implicitly on di, we finally 
obtain 

X{V3,V4) =P{9^)A{vs,V4)Q-\9^) 

^ (3.24LB) 
Y{vi,V2)=Q{ei)B{vi,V2)P-\ei) 

Here all the matrices depend on 61,62, 63 only via the arguments explicitly shown. 

Note that (3.23) imply the commutativity relation (3.4) for the modified model. (Using 
then simple arguments like those at the end of Sect. 2 of ref.[10] we can easily extend (3.4) 
to the original model as well). 

Thus, the relations (1.15)-(1.20) imply the star-star relation (2.36) which in its turn 
imply (3.20) and the relation (3.24). 

From (1.2), (3.3) and (3.24) it follows that 

Z=p^TriAiv3,V4)B{v^,V2)^ 

=P^Y.^A,,{y3,v,)B,,{v,,V2)r 
j 

where the j-summation is over all the diagonal elements Ajj and Bjj of A and B. Assuming 
the largest term in the summation is unique and writing a' for Aqq and U for Bqq (choosing 
J = for the largest term), it follows that from (3.25), (2.25) 

= pa{v3,V4)h{vi,V2), (3.25) 

where p is given by (2.23). This is the factorization property which will be used in the 
next Section. 
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3.2 Partition function per site. 



If L is infinite in all directions, we can evaluate the partition function per site solely 
from the symmetry and the factorization properties. In fact, the definition (1.1) is invariant 
under the symmetry transformation of the lattice generated by the symmetry group of the 
cube Q{C). On the other hand from (2.31) and (2.33) the effect of any such transformation 
on the weight function V is to multiply it by some equivalence transformation factors 
and permute cto, cti, a2, eta- Since two such transformations in (2.31), (2.33) generate all 
possible permutations of ciO: cii, ci2, 0:3, the partition function obey the following symmetry 
relation 

k(q;o, Q!i, q;2, 0:3) = K{a.i^ aj, ak, ai) (3.26) 

where /c, /) is any permutation of (0, 1, 2, 3). 

Thus we have two functional equations for k, the factorization property (3.25), (2.23) 
and the symmetry property (3.26). We shall see that these two functional equations define 
K up to multiplicative constant which is determined then from the normalization (2.29). 

The partition function per site k obviously depends on the number of the spin states 
N. We can exhibit this writing it as kn- Let us concentrate on the A/"-dependence in the 
functional equations. The symmetry relation (3.26) does not involve N at all. However, 
there are two possible sources of the A'^-dependence in the factorization equation (3.25). 
First p in (3.25) depends on as given by (2.23). Second, the formulae (2.13) expressing 
the variables f 1, . . . , f4 through the angles 9 1,62, 6*3 have the N-dependence coming from 
(2.11). This, however, can be readily removed since the A'"-th powers of the variables 
vi, . . .,V4 are the same for any N (the N-th. powers of the RHS's of (2.13) do not have any 
A'"-dependence) . So redefining the functions a and b we we can rewrite (2.25) in the form 

K ^ pa{x3, X4)b{xi, X2) (3.27) 

where Xi = , z = 1, . . . , 4, 

xi = e'""^ tan(^i/2) tan(^2/2), X2 = -e'""^ cot(6'i/2) tan(^2/2) 

(3.28) 

X3 = e-*"nan(6'i/2)tan(^2/2), a;4 = -e-'"^ cot(^i/2) tan(^2/2) 



are expressed merely through the angles 9 1,62, ^3 and the only A'"-dependence in the RHS 
of (3.27) comes from the factor p. Writing kjv in the form 

= A^(V^)2(^-^)/^ (3.29) 



where Aat is a numerical constant and substituting it into (3.27) we get the factorization 
equation for 

i> = (2^)a(x3, X4)b(xi, X2) (3.30) 

where ^ is given by (2.24). The symmetry relation for ip remains, clearly, the same as (3.26). 
Comparing (3.29) with eq.(7.33) of ref.[18] we see that ip satisfies the same factorization 
and symmetry relation as those for the partition function per site of the Zamolodchikov 



17 



model, (which is the N = 2 case of the considered model). It was shown in [18] that 
these relations determine K2 to within a numerical constant, so we can conclude that 



ip = \k2 (3.31) 

If we define the function G{f3) by 

1 /"^ r -K X ^ 

\ogG{(3) = — / xcota;+ -tan- -log(2sina;) da; (3.32) 
27r Jo L 2 2 J 

then the result of [18] can be written as 

«2 = ^G(/3o)G(/3i)G(/32)G(/33) (3.33) 

where /^O: A, /32, /^s are given by (2.11). Note that (3.33) satisfies the normalization (2.29). 
Taking this into account and substituting (3.31) into (3.29) and then into (2.29) we obtain 

log«^(^l,^2,^3) = ^^^^^log«2(^l,^2,^3) (3.34) 

which seems to be very simple and interesting result. 
4. Free Boson Model 

It is interesting to learn whether the method described above is unique to the considered 
model or it can be applied to the other models as well. In the other words how many 
solutions are there for the relations (1.15)-(1.20)? At the moment we do not know an 
answer to this question but can give one more example. It is the free boson (Gaussian) 
model on the three dimensional lattice. 

Let now the spins at the sites of the lattice take continuous real values — cxd < s < cxd 
and the summation over each spin s is replaced by the integration . . .ds. With these 
modifications the definitions (1.1), (1-12) remain valid, but the weights w and 7 in (1.12) 
are now given by 



-%(J?V 



w{v, a) = exp ( 2(^1 -y) ) ' " exp(-ia6), $(a) = exp(-ia /2). (4.1) 



Obviously, 

If we define v and A(t!) as 



w{v,Q)^l. (4.2) 



v^l/v, /\{v)^l-v, (4.3) 
for any v then the relations (1.15), (1.17), (1.18) remain valid while (1.16) is replaced with 

/oo 
7(a,6)7(-6,c)d6 = 27r5(a-c). (4.4) 
-00 
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Further, if we also replace the spin summation by the integration in (1.19)-(1.21) then 
these formulae become valid as well (provided one formally sets = 1 in (1.4), (1.6), i.e. 
A{v) is now given by (4.3) and u) — 1). The scalar functions 0i and 02 are now given by 



V2 - Vl 



1-1/2 



01(^1, ^"2) 



27rz(l — vi){l — V2) 



1-1/2 



(4.5) 



V3 - V4 



27rz(l — V3)(l — V4) 



Let us parameterize f 1, . . . , f4 through the angles 6*1, 6*2, 6*3 setting Vi — Xi, i — 1, . . . , 4, 
where a;'s are given by (3.28). Then, in particular, we have (cf. (2.22)) 



where ^ is given by (2.24). Finally the RHS of (2.28) is replaced by S{so — S2 — S2 + sq) 
and we still have (2.26). 

We have now all the required relations to reproduce the results(3.26), (3.27) for the free 
boson model and calculate its partition function per site . In fact all the calculations are 
identical to those we did before except that the exponents in (4.7) and (2.23) are different 
and the factors in (1.16) and (2.22) are replaced by 27r in (4.4), (4.6). (This factors 
are mutually cancel and therefore have no effect on the final result). Taking into account 
these minor differences we obtain 



A nice feature of the Gaussian model is that its partition function can be found by the 
elementary "Gaussian integration" directly from the definition (1.1). Together with (3.33) 
this means that the partition function per site of the A^-state model kn can be expressed 
as a rational power of a free boson determinant. We see no obvious reasons of this fact 
and it would be quite interesting to learn why it is so. We postpone the evaluation and 
the analysis of this determinant to the separate publication. 
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(4.6) 




(4.7) 



log Kb = -21ogK;2 



(4.8) 
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Figure Captions 



Fig.l. This figure explains the approach used in this paper. 

Fig 2. Arrangements of the spins sq, • • • , sq on the corner sites of an elementary cube 
of the simple cubic lattice C; do, di, d2, ds denote four oriented diagonals of the cube. 

Fig. 3. A typical elementary cube of £, with corner spins SO)---)^o and the center 
spin a. 
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